Giant spin-dependent photo-conductivity in GaAsN dilute nitride semiconductor 
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A theoretical and experimental study of the spin-dependent photoconductivity in dilute Nitride 
GaAsN is presented. The non linear transport model we develop here is based on the rate equations 
for electrons, holes, deep paramagnetic and non paramagnetic centers both under CW and pulsed 
optical excitation. Emphasis is given to the effect of the competition between paramagnetic centers 
and non paramagnetic centers which allows us to reproduce the measured characteristics of the spin- 
dependent recombination power dependence. Particular attention is paid to the role of an external 
magnetic field in Voigt geometry. The photoconductivity exhibits a Hanle-type curve whereas the 
spin polarization of electrons shows two superimposed Lorentzian curves with different widths, 
respectively related to the recombination of free and trapped electrons. The model is capable of 
reproducing qualitatively and quantitatively the most important features of photoluminescence and 
photocurrent experiments and is helpful in providing insight on the various mechanisms involved in 
the electron spin polarization and filtering in GaAsN semiconductors. 



I. INTRODUCTION 

The dependence of the recombination time on the rela- 
tive spin orientation of photogenerated carriers on para- 
magnetic centers, namely the spin dependent recombina- 
tion (SDR) has been known for over 50 years since the 
first optically detected magnetic resonance (ODMR) ex- 
periments carried out by Geschwind et al. 1 in AI2O3. 
It has been observed in surface centers on crystalline 
silicon 2 -! 3 - and on several systems including dislocated 
silicon 4 , amorphous silicon- and later on in (Al)GaAs^— . 
In dilute nitride GaAsN samples SDR has been evidenced 
with record high values at room temperature in photolu- 
minescence (PL) experiments^— and recently in pho- 
toconductivity (PC) measurements^ 4 -, stimulating new 
applications of the SDR as a light- or electron spin- 
polarization detector. Furthermore it has been shown by 
ODMR that it is in fact Nitrogen-induced Gallium self 
interstitial deep paramagnetic centers, and not Nitrogen 
itself, which harness the SDR in GaAsN alloys^ 5 .. 

The key point in the SDR effect is the existence of deep 
centers inside the gap possessing an unpaired electron at 
equilibrium (paramagnetic centers). As a consequence, 
if a photogenerated electron in the conduction band and 
the resident electron on the deep center have the same 
spin, the photogenerated electron cannot be captured by 
the center (it is generally assumed that the triplet levels 
are not bound). On the contrary, when the photogener- 
ated electrons and the the deep center resident electron 
have antiparallel spins, the capture is efficient. There- 
fore the recombination time of photocreated electrons de- 
pends on the relative orientation of the free electron and 
of the unpaired electron resident on the center. After 
the capture the centers annihilate one of the two elec- 
trons through a spin-independent recombination with an 
unpolarized hole. Under the influence of linearly polar- 
ized light an equal number of spin up and spin down 



conduction band electrons are created which are likewise 
captured by the paramagnetic centers. However, under 
circularly polarized light one given conduction band elec- 
tron spin orientation is preferred over the other. This 
not only creates spin polarized conduction band electrons 
through the usual selection rules but also dynamically 
polarizes centers. By both removing conduction band 
electrons with opposite spin orientation and by blocking 
the equally polarized ones, this mechanism gives rise to 
very high and very long living conduction band electron 
spin polarization and to larger PL intensities under cir- 
cularly polarized light compared to those under linearly 
polarized excit at ion 6 -! 7 - ' 1 1 1 1 . 

The theoretical work on SDR can be divided into two 
groups: stochastic Liouville equatio n 16 ! 17 , mostly de- 
voted to silicon, and rate equations 6 -! 7 - ! 10 ' 15 ! 18 ! 19 which 
have been used to study SDR in silicon as well as in 
GaAs(N). The former focuses mainly in the process of for- 
mation and dissociation of electron pairs in centers while 
the latter studies the conduction band carrier recombi- 
nation into centers and valence band states. Haberkorn 
et altH studied the process of formation and dissociation 
of paired spin-up and spin-down electrons via the solu- 
tion of the stochastic Liouville equation including the 
Zeeman part of the Hamiltonian and a stochastic linear 
term. Similar models also include interaction terms be- 
tween the electron spins. These models are suitable to 
the understanding of ODMR and similar experiments. 
Weisbuch et al& and Kalevich et al*&&^ studied the 
recombination dynamics of carriers through a rate equa- 
tions model that considers spin up and spin down elec- 
trons coupled to a distribution of paramagnetic centers 
constituting the SDR channel. In later approaches an 
equivalent set of rate differential equations for the spin 
variables of electrons and centers as well as the popu- 
lations of free and trapped carriers was propose d 12 ' 13 , 
lacking however the recombination through non param- 
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agnetic centers which is important to interpret the power 
dependence of SDR phenomena. 

In this paper we present a comprehensive theoret- 
ical and experimental work on SDR related phenom- 
ena in GaAsN, mainly focusing on photocurrent exper- 
iments. We have developed a new non-linear transport 
model based on rate equations for electrons, holes and 
deep paramagnetic centers in GaAsN based on the ini- 
tial model proposed to explain the strong SDR effects 
measured in photoluminescence experiment a 12 ' 13 ' 15 i 18 ' 20 . 
We show that the electron recombination to non- 
paramagnetic centers, which was neglected in the first 
models, is essential to understand the SDR dependence 
as a function of the excitation power. The proposed the- 
oretical model is able to reproduce accurately the most 
important features of both PC and PL experiments. Fi- 
nally we have measured and modeled the variation of 
the spin-dependent photoconductivity and photolumi- 
nescence signals as a function of an external magnetic 
field in Voigt configuration. As expected, the SDR ef- 
fect decreases when the amplitude of the external field 
increases. The spin-dependent photoconductivity signal 
exhibits clearly a Hanle like curve, well described by the 
theoretical model. 

The paper is organized as follows. The next section 
presents the samples' characteristics and the experimen- 
tal setup. We begin by discussing a previously formulated 
rate equations mode l 12 ! 13 in Sec. IIIII which is suitable to 
understand PL experiments. Starting from this system 
of differential equations we develop the model to include 
non SDR channels in the form of non paramagnetic cen- 
ters in Sec. 1111 Al and a uniform external magnetic field in 
Sec. 1111 Bl These preliminary parts put the fundamentals 
for the interpretation and modelization of the SDR medi- 
ated transport phenomena developed in Sec. IIII CI Here 
we put forward a nonlinear system of partial differential 
equations by introducing drift terms for the conduction 
band electrons and valence band holes. Sec. IIVI begins 
with a brief description of the numerical methods used to 
solve the system of partial differential equations obtained 
in the previous section. We discuss here two main results 
on transport experiments and calculations: the SDR as a 
function of the incident laser power and applied magnetic 
field. Our main results are summarized in Sec. [V] 



II. SAMPLES AND EXPERIMENTAL SET-UP 

The samples under study were grown under the same 
conditions by molecular beam epitaxy at T = 410°C on 
(001) semi-insulating GaAs substrates, after a 400nm 
GaAs buffer layer. They consist of GaAsN semicon- 
ductor epilayers of different Nitrogen content and thick- 
ness. Sample A has a nominally undoped 50nm thick 
GaAso.g8iNo.019 epilayer, sample B has a nominally un- 
doped lOOnm thick GaAs0.9924N0.0076 epilayer and sam- 
ple C a 50nm thick Si modulation doped (doping den- 
sity 10 18 cm' 3 ) GaAso.979No.021 epilayer. The growth 



was terminated with a lOnm GaAs cap layer and no post 
growth rapid thermal annealing was performed. We have 
observed similar effects in other doped or undoped sam- 
ples with N composition varying in the range 0.7% to 
2.6%. 

The excitation light was provided by a Ti:Sa laser, 
propagating along the direction normal to the sample 
surface, either in mode-locked regime, yielding the gener- 
ation of 1.5ps pulses at a repetition frequency of 80MHz, 
or in continuous wave operation, at a wavelength X exc = 
840nm. The laser was focused to a ~ 150/im or ~ 40//m 
diameter spots respectively in between two Ag electrodes 
deposited onto the sample surface ~ 0.8mm apart or be- 
tween two ohmic contacts of variable shape and dimen- 
sion defined by lithographic techniques and contacted by 
ultrasonic wire bonding. The laser light used for the 
photo excitation is either circularly polarized (right a + 
or left cr~) or linearly (a x or a Y ) polarized along a direc- 
tion perpendicular to the surface of the sample (k) and 
modulated by a mechanical chopper. The sample con- 
ductivity was then measured using a synchronous lock- 
in amplifier from the voltage drop at the terminals of a 
lOkH, load resistor placed in series with the sample. A 
constant voltage in the range < V < 12 volts was ap- 
plied between the sample electrodes. The PL intensity 
was simultaneously measured with the same lock-in tech- 
nique by recording its total intensity filtered of the laser 
scattered light, and integrated by a silicon p-i-n photodi- 
ode. All the experiments have been performed at room 
temperature. 



III. MODEL 

In GaAsN samples illuminated by linearly polarized 
light, identical populations of spin-up and spin-down 
electrons are photogenerated in the conduction band and 
are subsequently captured by paramagnetic centers con- 
taining an electron with opposite spin. The newly formed 
pairs dissociate at the centers living the same popula- 
tion of spin up and spin down paramagnetic centers as 
well. On the contrary, under circularly polarized excita- 
tion er~ approximately 75% of the photogenerated 
electrons are spin-up (spin-down) polarized and 25 % are 
spin-down (spin-up) polarized given the optical selection 
rules in GaAsSi. Thus, most of the processes occurring 
in the SDR channel will correspond, under a + excitation, 
to spin down conduction electrons rapidly recombining in 
spin up paramagnetic traps living most of the spin-down 
centers unpaired. After a few cycles the population of 
spin-up traps will be considerably depleted. This has the 
effect of a double spin filter: first, photogenerated spin- 
down electrons are prevented from recombining and sec- 
ond, electrons which undergo spin flipping mechanism are 
rapidly removed from the conduction band by the traps 
owing to the high availability of matching centers. More- 
over, the minority spin-up centers can efficiently trap a 
spin-down electron increasing the centers' spin polariza- 
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tion, so that the whole spin system is self-containing. 
At this stage the SDR channel is almost closed and the 
electrons-hole recombination through band to band op- 
tical transitions becomes more probable, increasing the 
photoluminescence yield. Moreover, as a result of the 
highly populated conduction band, a higher photocon- 
ductivity is expected under circular polarization com- 
pared to linear excitation. 

We start with the rate equations model originally pro- 
posed by Kalevich et a? i 12 ' 13 i 18 which grasps the essential 
features of the SDR effect: 
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Here the subindindices + and — stand for spin polarized 
electrons with spin projection +H/2 and —h/2 respec- 
tively along the incident light axis, thus n + and n_ are 
the number of spin-up and spin-down conduction band 
electrons of total number n = n + + n_. Similarly N + 
and N— are the the spin-up and spin-down paramagnetic 
traps. The total number of unpaired paramagnetic traps 
is given by Ni = N + + N_ and N 2 is the total number 
of electrons singlets hosted by the paramagnetic traps. 
Holes (p) are considered unpolarized as their spin relax- 
ation time is shorter than lps at room temperature^ due 
to the fourfold degeneracy of the valence band at k = 
and the large spin-orbit coupling in III-V semiconductors. 
Eqs. ©-(Ill) ensure conservation of charge neutrality and 
number of centers: 



n - p + N 2 = 0, 
Ni + N 2 = N. 



(5) 
(6) 

The terms of the form —"f e n±N^ are responsible for 
the spin dependent free electron capture in paramag- 
netic centers with recombination rate j e while the terms 
—jh,pN 2 model the spin independent recombination of 
one electron of the paramagnetic center singlet with a 
hole. The structure of these terms guarantees that spin- 
up conduction band electrons will only be trapped by 
spin-down centers and viceversa in order to comply with 
the Pauli exclusion principle. This recombination is at 
the heart of SDR mechanism: under cr + polarized exci- 
tation, for instance, the iV_ population builds up and 
the N + population is symmetrically quenched due to 
the spin dependent electron recombination on paramag- 
netic centers. Values of the spin relaxation time of free 
and unpaired electrons on the centers t s = 130ps and 
t sc = 2000ps respectively and the effective hole life time 
Th = 35ps as well as the typical ratio of the electron to 
the hole recombination coefficients ^e/lh = 10 ( where 
7/j = 1/rhN) are estimated from previous time resolved 
PL experiment s 10 ! 18 . 



The photo generation of electrons is accounted for by 
the terms G+ and G_ and of holes by G = G + + G_. Un- 
der CW excitation the generation term per unit surface 
is taken as 
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where P is the laser irradiance, hv is the photon's energy, 
2So is the maximum polarization degree of the generated 
electrons, a is the absorption coefficient of GaAs and L 
is the width of the GaAsN layer. The excitation light po- 
larization is gradually switched from linear (So = 0) to 
circular (So = 1/2) by the tanh function; the St parame- 
ter specifies the switching time between the the linearly 
(S z =0) and circularly polarized excitation (S z = So)- 
given the optical selection rules for GaAs under a non res- 
onant excitation involving only the valence band s = 3/2 
levels, the polarization of the generation term is set to 
So = 1/4 when the polarization of light is 100% circular. 
Under pulsed excitation the generation terms are given 
by 
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where P is the laser's time average irradiance for pulses 
of period t pp and width St- In our experiments t pp = 12ns 
and St = 1.5ps. 

In CW conditions the total photoluminescence inten- 
sity under linear (X) or circular (+) excitation is calcu- 
lated as Ix,+ = li-n (t) p (t) where 7 r is the optical recom- 
bination rate and t is a sufficiently long time to insure 
that the steady state has been reached. Under pulsed 
excitation the photoluminescence intensity is calculated 
as the time average Ix,+ = 7r f Q tpP n (t)p(t) dt/t pp . The 
total PL acts only as a probe of SDR effect here, since it 
weakly modifies the n and p densities^. 

A convenient pointer of SDR related phenomena is the 
SDR ratio (SDR r ) which is given by the quotient of the 
total photoluminescence intensity under circularly polar- 
ized excitation to the total photoluminescence intensity 
under linearly polarized excitation. The photolumines- 
cence SDR ratio is thus defined as 



SDR r = --. 
Ix 



(9) 



In the absence of the SDR mechanism, SDR r = 1 whereas 
in its presence SDR r > 1. In photocurrent experiments 
we replace the photoluminescence intensity by the pho- 
tocurrent in the definition of SDR r given that larger pho- 
tocurrents are expected in the presence of circularly po- 
larized light. 



A. Non-SDR channels: Non paramagnetic centers 

Although the main features of the SDR effect in di- 
lute Nitrides such as the very long conduction electrons 
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FIG. 1. (color on line) Photoluminescence SDR ratio vs laser 
power P for various ratios of the paramagnetic to non para- 
magnetic center of N/Ni — 4 (purple), 6 (blue), 8 (green), 
9 (red) and 10 (black). The (red) squares indicate the ex- 
perimental results. The dependence of the SDR maximum 
SDR m(la: as a function of N/Ni is shown in the inset. 



effective spin life time (in time resolved photolumines- 
cence experiments) 15 ' 18 are well described by the model 
above, it fails to reproduce the power dependence of the 
SDR ration. Fig. Q] shows the photoluminescence mea- 
surement [(red) squares] for sample A in the pulsed exci- 
tation regime and theoretical calculations [(black) solid 
line] of the SDR r as a function of the laser power. The 
data shows a steep increase of the SDR r as a function 
of the incident power, up to an optimum value. Above 
this power, the SDR r smoothly but monotonously de- 
creases. This is clearly due to the competition between 
SDR and non-SDR channels. In the low power regime 
(P < 50mW) the SDR mechanism dynamically polarizes 
the spin of paramagnetic centers thus blocking the re- 
combination of spin polarized conduction band electrons 
with the same spin orientation and filtering out oppo- 
sitely spin polarized ones. Once the optimum power is 
reached (P ~ 75mW in Fig. [1]) most of the paramagnetic 
centers are spin polarized. For higher laser powers the 
SDR mechanism is overwhelmed by the larger number of 
generated conduction band electrons and the recombina- 
tion through non paramagnetic centers dominates. 

This competition mechanism is a non linear contri- 
bution that has an important role at high excitation 
power and was not considered in the first models de- 
veloped to interpret the remarkable spin properties of 
GaAs N 12 ' 13 ' 18 . The presence of such centers is at- 
tested by the observation of numerous peaks in infrared 
spectroscopy^. Introducing this mechanism into the pre- 
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Here Ni and N 3 are the number of total and occupied 
non-paramagnetic centers respectively; 7 a is the recom- 
bination rate of conduction band electrons to non para- 
magnetic centers and 75 is the recombination rate of non- 
paramagnetic centers to valence band. 

Assuming that the non paramagnetic centers have the 
same recombination rates as paramagnetic centers (i.e. 
7<j/7e ~ 1 and 76/7/1 ~ 1), the main adjustable parame- 
ter of our model is the total number of non paramagnetic 
centers. Their role is clearly evidenced in Fig. [1] In this 
graph we present a series of SDR r vs power curves for var- 
ious rates of paramagnetic centers to non paramagnetic 
centers N/Ni keeping TV constant. As N/Ni increases, 
the SDR r maximum increases and moves to higher irradi- 
ance values. As the number of non paramagnetic centers 
is increased, the effective conduction electron recombi- 
nation time l/iaNi [see Eq. 1]10[) ] is shortened making 
this mechanism more likely at lower laser irradiance. At 
the same time the SDR maximum ratio is reduced as 
a consequence of the increasing transit of non polarized 
electrons through the non paramagnetic centers 15 . The 
best agreement with the experimental results for sample 
A (Fig. [TJ is achieved with N t = N/10 where the total 
paramagnetic center density N — N + + N- + N 2 is kept 
constant. The fitting of the PL SDR r power dependence 
allows us to fix the ratio N/Ni to be used in the following 
photoconductivity model. 



B. Effect of a magnetic field 

In order to include magnetic field effects into the 
model, we first recall the simple case of a spin 1/2 particle 
coupled to the magnetic field via the Zeeman Hamilto- 
nian 



H = fiu) ■ s 



(15) 



where us = g^B^B, g is the gyromagnetic factor of free 
electrons, [Ib is the Bohr magneton and s is the dimen- 
sionless vector spin operator for the particle. By solving 
the von Neumann equation for the density matrix and 
averaging the spin vector one finds that the Hamiltonian 
(1151) yields the well known coherent term contribution for 
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the total spin evolution 
dS 



dt 



u) x S. 



(16) 



The spin component along the propagation of light direc- 
tion is related to the number of spin-up and spin-down 
electrons in h units through S z = (n + — n_) / 2 12 i 13 . 

For the trapped electrons we can work out a similar 
dynamical equation given by 



dSc 



= fi x S r - 



(17) 



where f2 = <? c /iB?i-B with g c the gyromagnetic factor for 
trapped electrons and the spin component along the z 
axis is S zc = (N + — N-) /2. The gyromagnetic factors 
for free electrons and trapped electrons were set to g — 
0.5 and g c = 2 respectively-"^^^. 

Subtracting Eq. (|T0| for spin-down to spin-up elec- 
trons and adding the term (|16[) for spin precession we 
obtain 
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where we have used the identity 

n+N_ + n_N+ = ± (nNi - AS Z S ZC ) 



(18) 



(19) 



Adding Eq. (|10[) for spin-up and spin-down electrons 
considering isotropic spin relaxation time t s we obtain 
the rate equation for the total number of conduction band 
electrons. An analogous procedure is followed for the 
paramagnetic center variables. 

Finally the whole set of rate equations for the spin 
variables, total conduction band and trapped electrons is 
given by 

S + y (SN X - S c n) + ^S (Ni - N 3 ) 



+ — S + S x u = AG, 
S c + ^ (S c n - SATi) + — S c 
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Sxfl = 0, (21) 



(22) 
(23) 



Ni + f (nJVi - AS ■ S c ) - lh N 2P = 0, (24) 



N 2 - ^ (nNi - AS ■ S c ) + lh N 2P = 0, 
N 3 = la {N -N 3 )n- y b pN 3 



(25) 
(26) 



should be invariant under rotations, hence we have re- 
placed nNi - AS Z S ZC by nN x -AS- S c . By doing this the 
recombination rate of conduction electrons to paramag- 
netic centers is increased when S and S c are antiparallel 
whereas it vanishes when they are parallel as expected 
from the Pauli exclusion principle needed to form a sin- 
glet state. 



C. Photoconductivity 

Having dchncd the most important generation terms 
we turn to the central problem of this work: electron 
transport combined with SDR effects. 

Modulation of the photoconductivity with the polar- 
ization of the excitation light up to 40% has been indeed 
measured very recently at T — 300 Kii. 

Several theoretical contributions have been presented 
to study the transport of electrons in doped GaAsN^l. 
Fahy et al. used a tight binding model and the Boltz- 
mann equation in the relaxation-time approximation to 
obtain the mobility of conduction band electrons in 
GaAsN^i. Nevertheless the spin degree of freedom was 
not considered in this work. 

In order to focus on SDR associated phenomena we 
present a simpler approach in which the electron and 
hole density currents are introduced through drift terms 
for the current neglecting other transport phenomena as 
diffusion. The Einstein-Smoluchowski relation for the 
diffusion coefficient yields a diffusion length of a few hun- 
dred nanometers compared to the laser spot size in the 
order of the 100/xm due to the low mobility (^i « 100 
cm 2 V/s — ~— ) of electrons and holes in GaAsN. There- 
fore diffusion currents are neglected in the calculation. 

In this case n± and p are understood as the ID local 
densities of electrons and holes, therefore n± = n± (t, x) 
and p = p(t,x). Similarly the number of paramagnetic 
impurities can be replaced by the local density of impu- 
rities N± = N± (t,x) and N 2 = N 2 (t,x). S = S(t,x) 
and S c = S c (i, x) are the free electron and paramagnetic 
trap spin densities respectively. 

The drift density current for electrons is J n = J + + J_ 
where J± = —fi e En± corresponds to spin up- and spin 
down-electrons. For holes the drift density current is 
Jp = HhEp. Here [i e and \ih are the electron and hole 
mobilities respectively. They are set to fi e — 200 cm 2 /Vs 
and Hh — 100 cm 2 /Vs respectively as measured and cal- 
culated for similar samples 2 ^—. We must also include 
the spin density current terms. According to Dyakonov 
et o/j^i, the spin polarization current density tensor, ne- 
glecting a small spin orbit contribution, can be written 
for electrons as 



(27) 



where AG = (G + — G_)fc/2. In an isotropic media 
the dynamical equations for the three spin components 



where D is the diffusion coefficient. Due to the small 
conduction electron mobility the diffusion term can be 
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neglected in Eq. 



27jl so that 

Js ~ ~HeES. 



(28) 



Gathering all the terms we obtain the following system 
of equations 
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(29) 
(30) 

(31) 
(32) 

(33) 
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The density of photogenerated carriers is obtained by 
multiplying (J7J or © by the Gaussian spatial distri- 
bution exp — (x — xq) 2 /2Sx /V2tts x which specify the 

laser profile needed in transport simulations. The laser 
spot radius s x is varied according to the experimental 
conditions. Even though the generation terms for trans- 
port carriers do not figure explicitly in the equations 
above, they are considered through the boundary con- 
ditions. The number of incoming electrons and holes on 
opposite sides of the sample are fixed for a given voltage 
and electrical current. 

It is worthwhile mentioning that, while the expression 
for the local conservation of paramagnetic centers [Eq. 
(|6])] remains unchanged, the charge conservation condi- 
tion [Eq. ([5])] is replaced now by the continuity relation 



(n-p + N 2 + N 3 ) + — (J e 



8t- " —-■Ox- ^ ) = °- (36) 
As for the PL experiments it is convenient to define 
a parameter that allows us to pin the SDR phenomena. 
The transport SDR ratio, similarly to the photolumines- 
cence SDR ratio, is thus the quotient of the time averaged 
density current in circular polarization (J+) to the time 
averaged density current in linear polarization (Jx)' 



SDR r = J±- 
Jx 



(37) 



where explicitly J+,x = f PP J+,x (t)dt/t pp . 

IV. RESULTS AND DISCUSSION 

The rate Eqs. (|2TH) - (|26I) were solved by fourth-order 
Runge-Kutta method. As initial conditions the amount 
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FIG. 2. (color online) (a) Photoluminescence intensity and 
(b) photocurrent signal as a function of time under circularly 
(<r + ) and linearly {a ) polarized light. The (red) squares 
indicate the experimental results for PL and PC and the thick 
black line represents the theoretical calculations. 



of spin-up and spin-down electron on paramagnetic cen- 
ters was set to N/2 while the remaining variables were 
set to zero. 

In the transport case the carrier generation terms were 
introduced through the boundary conditions for a fixed 
voltage and current intensity; Dirichlet boundary con- 
ditions were imposed separately at the two ends of the 
sample for electrons and holes given that their distribu- 
tions build up at opposite extremes. Accordingly, the 
space was discretized by forward and backwards differ- 
ences for electrons and holes respectively. Initially the 
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FIG. 3. (color online) (a) Spin-up n+ (green), spin-down n_ 
(red) and total n (blue) density of free electrons, (b) Density 
of trapped electrons on paramagnetic centers: Spin-up N + 
(green), spin-down N- (red) and Ni = N+ + N- (blue). 



spin S c of trapped electrons was set to zero uniformly 
along the sample. As the mobility in GaAsN samples 
is very low it is possible to neglect the effects of trans- 
port in a first round of calculations using fourth-order 
Runge-Kutta method for each discretized position. The 
resulting system was solved by point Gauss-Seidel itera- 
tion method including the transport terms. In this cal- 
culations we used the same fitting parameters as the one 
introduced in Sec. IIIII 

In Fig. [5] we present the simultaneous measurement 
of photoluminescence (a) and photocurrent (b) in sam- 
ple B with lithographic contacts under CW excitation of 
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FIG. 4. (color on line) Photocurrent SDR r vs laser irradiance 
P for different values of the magnetic field. The (red) squares 
are the experimental results, the dark solid lines indicate the 
theoretical curves for different transverse magnetic fields OmT 
(black), 200mT (red) and 400mT (green). 



P = 5mW. The (red) squares indicate the experimen- 
tal measurement of PL and photocurrent while the solid 
lines correspond to the theoretical calculations obtained 
from the solution of the set of differential Eqs. (l^ - d^SI) . 
In this graph we observe an increase in the PL intensity 
under circularly polarized light compared to a lin- 
early polarized excitation {c x ). Likewise the photocur- 
rent increases under a a + excitation with respect to the 
a x case. A clear correlation is observed between the PL 
and PC behaviors. Under a x excitation the density of 
spin-up and spin-down polarized electrons are identical, 
namely n + — n_ whereas with a a + excitation n_ be- 
comes considerably larger than n + [Fig. [3] (a)]. Through 
the dynamical polarization mechanisms the iV_ to N+ 
population ratio correspondingly follows the conduction 
band electron one [Fig. [U(b)]. When ?i_ ^> n + , the JV_ 
density sharply increases, blocks further recombination of 
conduction band electrons of same spin orientation, and 
the iV_|_ population is depleted. Accordingly, for linearly 
polarized light we can see identical N + and A_ popula- 
tions. The same model reproduces well the dependence 
of the PL and PC signal as functions of the excitation 
light polarization [solid line in Figs. [U (a) and[D(b)]. 

Fig. Q] [(red) solid squares] shows the photocurrent 
SDR r as a function of power P under pulsed laser exci- 
tation for sample C whereas the solid (black) line repre- 
sents the theoretical photocurrent SDR r obtained from 
the set of differential Eqs. ([!?!))) - (|55j) . Features identical 
to the PL experiments are measured and are precisely 
reproduced by the model. 

In order to get further insight we studied the Hanle 
effect i.e. the depolarization of electrons in the presence 
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FIG. 5. (color on line) Hanle curve, (a) Photocurrent SDR 
ratio vs magnetic field B. The (red) squares display the exper- 
imental results for a power of 80mW in pulsed excitation and 
the (red) circles exhibit the experimental results at lOOmW 
in CW excitation. The solid lines correspond to the theo- 
retical calculations, (b) Theoretical values of the conduction 
band electron circular polarization P c vs magnetic field B for 
P = 80 mW (black), 40 mW (red) and 20 mW (green) in 
pulsed operation. 



of a magnetic field B in Voigt configuration^ and the 
consequent quenching of the SDR rate. For non vanish- 
ing values of the magnetic field (red for 200mT and green 
line for 400mT) a progressive decrease in the SDR r is ob- 
served in Fig. @] as expected. As the transverse magnetic 
field is increased the originally polarized spin of conduc- 
tion band and trapped electrons begin to precess around 
the magnetic field with different angular velocities given 
the different g factors. The recombination through para- 
magnetic centers is thus triggered due to the dephasing 



induced mismatch between spin of conduction band elec- 
trons and trapped electrons. 

In Fig. [5] (a) we present the experimentally recorded 
variation of the SDR ratio as a function of an external 
magnetic field for sample C at P — 80 mW under pulsed 
excitation conditions. A very strong reduction of the 
photocurrent SDR ratio describing a quasi Lorentzian 
curve is observed. The model described in Sec. MI Bl 
reproduces accurately the measured magnetic field de- 
pendence [(black) solid line in Fig. [5] (a)]. To clarify 
further the magnetic field dependence let us recall that 
in the simplest situation of an electron spin precessing 
around the external magnetic field, the Hanle effect yields 
a Lorentzian dependence for the electron spin polariza- 
tion described by: 



Pr.= 



where the half- width -B1/2 is given b y 21 i 32 

Bl/2 = 



(38) 



where T s is the mean decay time: 



1 



1 

Ts 



1 
T 



(39) 



(40) 



and t is the electron lifetime. We emphasize here that 
the data displayed in Fig. [5] (a) do not correspond to 
a classical Hanle curve which probes the spin depolar- 
ization in a photoluminescence experiment. Therefore, 
to interpret the experimental data, we have calculated 
in Fig. [5] (b) the free electron spin polarization degree 
as a function of the magnetic field B obtained from the 
model for various values of the power for circularly po- 
larized excitation. This curve presents new features 
hidden in the SDR magnetic field dependence in PC. The 
curves show two superimposed Lorentzian with different 
half widths that strongly suggest two distinct processes 
with different mean spin life times. The magnetic field 
has here the double effect of depolarizing the conduction 
band and paramagnetic center resident electrons. Due to 
their significantly different spin lifetimes the B field de- 
pendence of the polarization shows the double Lorentzian 
curve as observed by Kalevich et al£2- in CW PL polar- 
ization experiments. For low excitation power the narrow 
Lorentzian is controlled by the localized electrons (long) 
spin life time while the larger Lorentzian corresponds to 
the free electron (short) spin life time. In the photocon- 
ductivity experimental results of Fig. [5] (a), the width of 
the curve is determined by the localized electron spin life- 
time T s which corresponds to the electron lifetime r. The 
latter is controlled by the spin-dependent capture of free 
electrons (see Sec. IIII|) . To further support this interpre- 
tation we have measured the magnetic field dependence 
of the photocurrent at much lower photo-generated elec- 
tron density (i.e. under CW excitation). As shown in 



9 



Fig. [5] (a) [(red) circles], the width of the curve is one 
order of magnitude smaller due to the longer localized 
electron lifetime^. We emphasize that the same set of 
parameters were used in the calculations for Fig. [3J Fig. 
Hand Fig. 

V. CONCLUSIONS 

We have presented photoluminescence and photocur- 
rent measurements in GaAsN samples that evidence the 
presence of many SDR related phenomena such as a spin 
dependent photoconductivity change and Hanle effect in 
the SDR ratio of the photoconductivity. We have found a 
strong correlation between photoluminescence and pho- 
toconductivity results. To get further insight on the in- 
terplay of the many different mechanisms involved we 
developed a model that correctly reproduces the main 
features of these phenomena. The model is based on 
a set of non-linear transport differential equations con- 
taining drift density currents for electrons and holes and 
terms that account for photogeneration, SDR in para- 
magnetic centers, and recombination through non para- 
magnetic centers. The presence of a constant magnetic 
field is introduced through spin-precession terms. First, 
we have shown that the interplay of SDR and non SDR 
channels is responsible for the most important features 



of the SDR ratio vs laser irradiance curve. It presents a 
maximum that can be explained as a result of the com- 
petition between recombination through SDR and non 
SDR channels that activate at high laser pumping. As 
the recombination through non paramagnetic centers be- 
comes more likely the maximum SDR ratio decreases. 
Second, we investigated the effect of a transverse mag- 
netic field (Voigt configuration) on the Spin Dependent 
Photoconductivity. The SDR ratio vs magnetic field 
curve shows a Lorentzian-like shape which is reminiscent 
of the Hanle effect experienced by conduction band elec- 
trons and traps. 
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